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SUSY on the lattice 

I. Montvay 

'^Deutsches Elektronen Synchrotron DESY, 
Notkestr. 85, D-22603 Hamburg, Germany 

The motivation and perspectives of numerical simulations of supersymmetric Yang-Mills theories are reviewed. 
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1. INTRODUCTION 

Supersymmetric quantum field theories, re- 
viewed for instance in refs. have peculiar 
properties. They represent an interesting special 
subset of quantum field theories. Because of their 
highly symmetric nature, they are best suited for 
analytical studies, which sometimes lead to exact 
solutions Q. 

1.1. Supersymmetric Yang-Mills theory 

Since local gauge symmetries play a very im- 
portant role in Nature, there is a particular inter- 
est in supersymmetric gauge theories. The sim- 
plest examples are supersymmetric Yang-Mills 
(SYM) theories which are supersymmetric ex- 
tensions of pure gauge theories. Depending on 
the number N of supersymmetry generator pairs 
QiaiQia, {i = 1,2,...,A^), One speaks about 

= 1 "simple" supersymmetry or iV > 2 "ex- 
tended" supersymmetry. Interesting examples of 
SYM theories have, for instance, TV = 1 or 2 su- 
persymmetry. 

The SYM action with = 1 supersymmetry 
can be written as follows: 

SsYM ^iinlyT j d'^x(feTT{w°'Wo:)^ 
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Here r is a complex coupling parameter 

^=2^ r 
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combining the gauge coupling g with the 0- 
parameter. The first line of eq. ([|) is in superfield 
notation, in terms of the spinorial field strength 
superfield W{x, 9, 6)a which depends on the four- 
coordinate X and the anticommuting Weyl-spinor 
variables Oa,Oa {oi,a — 1,2). After performing 
the grassmanian integration on 9, one obtains the 
second form of the SYM action in eq. (|l|) in terms 
of the component fields: the field strength tensor 
Ffj_^, its dual F^^ = \e^vpaFpcr, the Majorana 
fermion field in the adjoint representation A, A 
and the auxiliary field D. Of course, denotes 
gauge covariant derivative and = {1, cr^'^''^}. 
(For details of the formalism see refs. 

Since the auxiliary field without kinetic term 
D can be trivially integrated out, Ssym is noth- 
ing else but a Yang-Mills theory with a massless 
Majorana fermion in the adjoint representation. 
Hence the Yang-Mills theory of such a field, called 
"gaugino" or, in the context of strong interac- 
tions, "gluino" , is automatically supersymmetric. 
Introducing a non-zero gaugino mass rrig breaks 
supersymmetry "softly" . Such a mass term is: 



(3) 



Here in the first form the Majorana- Weyl compo- 
nents A, A are used, in the second form the Dirac- 
Maj orana field ^. 

Ssym in eq. defines the simplest, = 1 
SYM theory. It is reasonable to start the non- 
perturbative lattice studies by this. An impor- 
tant feature of this theory is its global U{l)x chi- 
ral symmetry which coincides with the so called 
R-symmetry generated by the transformations 
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Q' — JVC 
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This is equivalent to 

A'„ = e'^^K , = e~'^\a , = e-^'^T=* . (5) 

The C/(l)A-synimetry is anomalous: for the cor- 
responding axial current = 4* 7;^ 75 4*, in case of 
SU{Nc) gauge group with coupling g, we have 

However, the anomaly leaves a Z2n^ subgroup of 
U{1)\ unbroken. This can be seen, for instance, 
by noting that the transformations 

^^6"*'^^=^*, ^^^le-'-f^^ (7) 

are equivalent to 

mg ^ m~ge-^'^^^ , Qsym ^ Qsym - 2N,ip , (8) 

where Osym is the 0-parameter of gauge dynam- 
ics. Since Qsym is periodic with period 27r, for 
nig — the U{1)\ symmetry is unbroken if 

<P = ¥'fc = ^, (A; = 0,l,...,2iV,-l) . (9) 

For this statement it is essential that the topolog- 
ical charge is integer. 

The discrete global chiral symmetry Z2n^ is ex- 
pected to be spontaneously broken by the non- 
zero gaugino condensate {XX) 7^ to Z2 defined 
by {vqtVnA (-^ ^ is ^ rotation). Instan- 
ton calculations (for a review and references see 
give at SsYM = in A^c degenerate vacua 
(fc = 0, . . . , iV, - 1) 

Here Asy m denotes the A-parameter of the gauge 
coupling. The transformation with Qsym is 

(AA)e,,„, =e^«— /^^(AA}o. (11) 

In general, an important and interesting ques- 
tion is whether supersymmetry can be broken 
spontaneously or not. In pure SYM theory, with- 
out additional matter supermultiplets, this is not 
expected to occur. An argument for this is given 
by the non-zero value of the Witten index j^] 

— Tl ( 1) n,}joson '^fermion : (1^) 



which is equal to the difference of the number 
of bosonic minus fermionic states with zero en- 
ergy. It is supposed not to change with the pa- 
rameters of the theory. For SYM theory we have 
wsYM — Nc, therefore there is no spontaneous 
supersymmetry breaking. (The Nc ground states 
discussed above correspond to this Witten index.) 

1.2. Seiberg- Witten solution 

The SYM theory with N = 2 extended super- 
symmetry is a highly constrained theory which 
has, however, more structure than the relatively 
simple = 1 case discussed above. In particular, 
besides the iV = 1 "vector superfield" containing 
the gauge boson and gaugino (A^,A), it also in- 
volves an = 1 "chiral superfield" {(f>, A') in the 
adjoint representation which consists of the com- 
plex scalar <j> and the Majorana fermion A'. The 
Majorana pair (A, A') can be combined to a Dirac- 
fermion ijj and then the vector-like (non-chiral) 
nature of this theory can be made explicit. 

The main new feature of A*" = 2 SYM is that 
it contains also scalar fields, hence there is the 
possibility of Higgs mechanism. Let us here 
consider only an SU(2) gauge group, when the 
Higgs mechanism implies the symmetry break- 
ing SU{2) U{1). The complex expectation 
value of the scalar field {(p) parametrizes the mod- 
uli space of zero-energy degenerate vacua. The 
degeneracy is a consequence of extended super- 
symmetry. Due to the Higgs mechanism the 
"charged" gauge bosons become heavy and the 
low-energy effective theory is A^ = 2 SYM with 
U{1) gauge group. 

As it has been shown by Seiberg and Wit- 
ten Q , extended supersymmetry and asymptotic 
freedom can be used to determine exactly the low- 
energy effective action in the Higgs phase, if the 
vacuum expectation value is large. It turns out 
that, besides the singularity which corresponds 
to asymptotic freedom at infinity, at strong cou- 
plings there are two singularities of the effective 
action which describe light monopoles and dyons, 
respectively. 

The unprecedented achievement of an exact so- 
lution for a non-trivial four-dimensional quan- 
tum field theory induced a considerable activ- 
ity among theoreticians. Recently there is a new 
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way of considering non-perturbative problems in 
quantum field theory. It is an interesting question 
how much the more conventional approach based 
on lattice regularization can contribute to these 
new insights. 

2. LATTICE FORMULATION 

The simplest supersymmetric gauge theory is 
TV = I SU{2) SYM. Massive gluinos (Majorana 
fermions in the adjoint representation) break su- 
persymmetry softly. As discussed in section |lH 
the supersymmetric point is at rrig = 0. 

2.1. Actions 

The Curci- Veneziano action of iV = 1 SYM 
is based on Wilson fermions. The effective gauge 
action obtained after integrating out the gluino 
field is given by 

Scv^PY, (l-\TrUp^-]-\ogdetQ[U] ,(13) 
pi ^ ' 

where the fermion matrix is 

and the gauge link in the adjoint representation 
is defined as 



(15) 



The Majorana nature of the gaugino is rep- 
resented in eq. (^3|) by the factor ^ in front of 
logdetQ. There is no problem with the sign of 
det Q because we have det Q = det Q, where 



g = 75g - 



(16) 



is the hermitean fermion matrix. From the rela- 
tions 



cgc^i = g^ , BQB-^ = g^ 



(17) 



with the charge conjugation matrix C and B = 
C75, it follows that every eigenvalue of Q and Q 
is (at least) doubly degenerate. Therefore, with 
the (real) eigenvalues Xi of Q, we have 



detg = J|A^ > 0. 



The doubling of the eigenvalues of Q and/or Q 
implies that for fermions in the adjoint represen- 
tation zero eigenvalues always occur in pairs, as in 
the continuum (This is different for fermions 
in the fundamental representation: see for in- 
stance Q.) 

Of course, the lattice action in eq. (|3|) is not 
unique. Another possibility is based on five- 
dimensional domain walls pO[ . In this approach 
one knows the value of the bare fermion mass 
where the supersymmetric continuum limit is 
best approached and one may have advantages 
from the point of view of the speed of symme- 
try restoration, but one has to pay with the pro- 
liferation of (auxiliary) fermion flavours. There 
are also speculations about direct dimensional re- 
duction on the lattice from c? = 10 or d = 6 di- 
mensions down to the physically interesting case 
d = 4, started by ref. ||ll[] and continued by 
ref. ||l2|. However, up to now it has not yet 
been demonstrated that this approach docs really 
work. 

2.2. Algorithms 

As one can see in the lattice action in eq. (p^, 
Majorana fermions effectively mean a flavour 
number Nf = ^. This can be achieved by the 
hybrid classical dynamics algorithm [ p^ , which is 
applicable to any number of flavours. Since it is a 
finite step size algorithm, an extrapolation to zero 
step size is needed. This approach has been im- 
plemented and successfully tested for the Curci- 
Veneziano action ( p^ ) by Donini and Guagnelli in 
ref. 0. 

Another possibility to simulate non-even num- 
bers of flavours is based on the multi-hosonic algo- 
rithm [ p5| . A two-step variant using a noisy cor- 
rection step ijl^ has been developed in ref. [p7| . 
This algorithm is based, for general Nf, on the 
approximation 

where the polynomial P„ satisfies 



lim P^{x) = x^^f'^ 



(20) 



(18) 



within a suitably chosen positive interval x G 
[e. A], which covers the spectrum of Q^Q in typi- 
cal gauge configurations. In the two-step variant 
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one uses a first polynomial Pn for a crude approx- 
imation and realizes a fine correction by another 
polynomial P„, where now 



approximated: x 



lim Pn{x)Pn{x) 



for X e [e, A] . (21) 



The fermion determinant is approximated as 



det(QTQ) 



1 



detPs(QtQ) detF„(gtQ) 



.(22) 



The advantage of this two-step variant is its 
smaller storage requirement and shorter autocor- 
relations, because the order of the first polyno- 
mial n can be kept small. (In praxis, for instance, 
n < 10—12 for lattice sizes < 8"^ • 16.) As a conse- 
quence of the smallness of n, this remains a large 
step size algorithm, also in the limit of good pre- 
cision. 

In order to achieve an exact algorithm in this 
two-step variant of the multi-bosonic algorithm, 
in principle, an extrapolation to n — > oo is 
needed, together with measurement-corrections 
for eigenvalues outside [e. A]. In praxis, however, 
{n,n) and [e. A] can be chosen such that these 
kinds of corrections are smaller than the statisti- 
cal errors. 

Random matrix models suggest that the fluc- 
tuations of the minimal eigenvalue are inversely 
proportional to the lattice volume. (For refer- 
ences and a recent summary see ref . pq ] . ) This is 
advantageous for the choice of the interval [e. A]. 
For the calculation of the necessary polynomials 
procedures written in Maple are available [l9| . 
The quality of the polynomial approximations is 
illustrated by figure |[ 

Besides the above "hermitean" version "com- 
plex" (non-hermitean) variants of the multi-boso- 
nic algorithms ||2^ can also be exploited, prefer- 
ably also in two-step variants. In order to im- 
prove performance, preconditioning according to 
ref. turned out to be very useful. 

2.3. Quenched computations 

Quenched calculations without the effect of dy- 
namical gluinos can be useful for roughly localiz- 
ing the critical region in the hopping parameter K 
and for testing the methods to determine masses. 
Such calculations have been performed recently 
in refs. |22| and p^. An example of the obtained 
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Figure 1. The deviation norm of the polyno- 
mial approximations of x^^^^ as function of the 
order for different values of A/e. The asterisks 
show the e/A hmit. (Ref. ||l|l) 

lowest masses is shown in figure ^. Notation con- 
ventions are as follows: tt and a denote, respec- 
tively, the pion-like and sigma-like states made 
out of gluinos. X is a Majorana fermion state 
made out of gluons and gluinos. 

3. FIRST RESULTS AND OUTLOOK 

The goal of the DESY-Miinster-Athens collab- 
oration is to study non-perturbative proper- 
ties of SYM theories by Monte Carlo simulations. 
The computations are performed on the APE- 
Quadrics at DESY-Zeuthen and on the CRAY- 
T3E at HLRZ-Jiilich by using the two-step multi- 
bosonic fermion algorithm for gluinos discussed in 
section 2.2. 



3.1. Chiral symmetry breaking 

An important first step in the numerical inves- 
tigations is to find the critical hopping parame- 
ter Kcr corresponding t o ze ro gluino mass. As 
it is discussed in section 



1.1 



the supersymmetric 
point is at zero gluino mass. This tells how the 
supersymmetric continuum limit f3 oo has to 
be performed in the {/3, K)--plane. The basic ex- 
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3 = 2.6 ; V = 16-^ • 32 
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Figure 2. Masses of tt, ct, x a-nd 2/9 as a function 
of 1/K (the mass of tt is squared). The extrapo- 
lation at Kcr is also shown (curves). (Ref. [ p3| ) 

pectation is that at K = Kcr, for gauge group 
SU{Nc), there are Nc degenerate ground states 
with different values of the gluino condensate (see 
eq. (|l0[)). This corresponds to the expected spon- 
taneous global chiral symmetry breaking pattern 

The coexistence of several different vacua im- 
plies a first order phase transition a.t K = Kcr, 
at least in the continuum limit. At finite lattice 
spacing in the {(3, if )-plane this might be repre- 
sented just by cross-over lines. Nevertheless, the 
simplest possibility is a first order phase transi- 
tion already for finite /3. For instance, in case of 
SU{2) gauge group the phase structure may be 
as shown by figure ||. An alternative phase struc- 
ture can be motivated by the modified Veneziano- 
Yankielowicz effective low energy action proposed 
recently by Kovner and Shifman ||2^. They sug- 
gest the existence of an additional massless phase 
with no chiral symmetry breaking. On the lat- 
tice this could lead to a more complicated struc- 
ture, for instance, because the symmetric phase 
could be stable or metastable in some intermedi- 
ate range of the hopping parameter. 

An interesting point is the dependence of the 
phase structure on the gauge group. For example, 
in case of SU{3) there are at least three degen- 
erate vacua and beyond the critical hopping pa- 
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Figure 3. Expected phase structure in the (/3, K)- 
plane. The dashed-dotted line is a first order 
phase transition at zero gluino mass. 

rameter K > Kcr we expect to be at Qsym = tt. 

Since the Monte Carlo simulations are usually 
done at non-zero gluino mass nig ^ 0, the best 
procedure is to derive analytically the dependence 
of different quantities on rrig ||2^ and compare 
them directly with numerical data. In case of the 
gluino condensate the best possibility is to study 
the difference A(AA) between the values of (AA) 
in different vacua. The advantage is that the ad- 
ditive ultraviolet divergences, due to the breaking 
of supersymmetry at finite lattice spacing, cancel. 

3.2. Bound state spectrum 

The low energy effective actions predict that at 
the supersymmetric point the low lying states are 
organized in massive supermultiplets. The sim- 
plest possibility is to have a lowest supermultiplet 
generated by a superfield made out of gluinos . 
On the lattice such states are tt, a and x discussed 
already in section 2.3. Other obvious candidates 
are glueball states observed in pure Yang-Mills 
theory, which might also belong to another super- 
multiplet. Of course, the notion of constituents 
is in general not expected to work perfectly in a 
strongly interacting theory. Therefore one has to 
be open minded towards non-trivial mixings. 

One of the main goals of the DESY-Miinster- 
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Athens collaboration is to determine the masses 
of low lying states in the vicinity of the supersym- 
metric point at zero gluino mass. First results are 
presented at this conference. In general, as shown 
by figure ^, the supersymmetric degeneracy of low 
lying states is not yet observed. Further simula- 
tions closer to Kcr and at larger f3 are in progress. 

Lattice: 6^x12, ^ = 2 J 
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Figure 4. Dependence of the lowest bound state 
masses on the hopping parameter K on 6^ ■ 12 
lattices at /3 2.3. Masses of states made out of 
gluons (g) and gluinos (g) are shown. (Ref. p^) 

3.3. Other questions 

Besides the phase structure and the low lying 
mass spectrum there are also other interesting 
questions which can be studied by Monte Carlo 
simulations. For instance, the direct investigation 
of broken supersymmetry Ward-Takahashi iden- 
tities is possible. In general, the restoration of 
supersymmetry in the continuum limit has to be 
understood. Working with gluinos will teach us 
a lot about the simulation of non-cvcn numbers 
of flavours and about fermions in the adjoint rep- 
resentation. We shall hopefully get answers to 
these and other related questions at forthcoming 
lattice conferences. 
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